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Abstract 



Using e + e~ annihilation data collected by the CLEO II detector at CESR, 
we have observed the decay Df — ► lutt + . Lhis final state may be produced 
through the annihilation decay of the or through final state interactions. 
We find a branching ratio of L(L>+ -> wvr + )/L(L>+ -> r?vr + ) = 0.16 ± 0.04 ± 
0.03, where the first error is statistical and the second is systematic. 

PACS numbers: 13.25.Ft, 14.40.Lb 
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It has been suggested that the uir + decay mode could be a clean signature for the 
annihilation decay of the 0. While the simple annihilation diagram can produce p°7r + , 
it cannot produce uir + , because this final state has isospin and G-parity I G = 1 + ; to do so 
would require a second-class axial current @. If at least two gluons connect the initial state 
quarks to the final state quarks, the decay — > ujtt + through the annihilation diagram 
is allowed. The possibility that this final state might arise through final state interactions 
(FSI) has also been extensively discussed Fermilab E691 set a 90% C.L. upper limit 

of F(D+ -»■ iV7t + )/T(D+ -»■ 07T+) < 0.5 0, or B(D+ -> uti+) < 1.8% 0; this is the most 
sensitive published limit. To date, the only clear evidence for the annihilation decay of a 
charmed meson is — > fi + u ||. This letter describes the first observation of the decay 
— > um + , and the measurement of the branching ratio T(D+ — ► cj7r + )/T(.D+ — > i]tt + ). 

A recent paper by Buccella et al. predicts nonresonant FSI should produce B(D+ — > 
D7Y + ) = 2.9 x 10~ 3 ||; however, their prediction for the related decay mode, Df — > r]'p + , 
does not agree well with measurements There could be a small contribution to the 

U7Y + decay rate from spectator decay, due to the tiny ss component of the uj. The ss content 
of the ui is estimated to be ~ 0.4%, assuming a vector octet-singlet mixing angle of 39° [[7]. 
The branching fraction for spectator decay to um + can naively be estimated to be about 
0.004 x B(Df — > </>7r + ) w 1.5 x 10~ 4 . This is below our current sensitivity. There may also 
be mixing of the u with the <fi through their common decay modes. 

The data used in this analysis were collected with the CLEO II detector |10| at the Cornell 
Electron Storage Ring (CESR). The detector consists of a charged particle tracking system 
surrounded by an electromagnetic calorimeter. The inner detector resides in a solenoidal 
magnet, the coil of which is surrounded by iron flux return instrumented with muon counters. 
Charged particle identification is provided by specific ionization (dE/dx) measurements in 
the main drift chamber. The data were taken at center-of-mass energies equal to the mass 
of the T(45) (10.58 GeV) and in the continuum approximately 50 MeV below the T(45). 
The total integrated luminosity was 4.7 fb _1 . 

Events used in this analysis were required to have a minimum of three charged tracks, 
and energy in the calorimeter greater than 15% of the center-of-mass energy. Charged tracks 
were required to have dE/dx measurements within 2.5 standard deviations of that expected 
for pions. Only energy clusters in the calorimeter with | cos#| < 0.71 (where 9 is the polar 
angle with respect to the beam axis) that were not matched to charged tracks were used as 
photons. Photons with energy greater than 30 MeV were combined in pairs to reconstruct 
7r°'s. The invariant mass of the two photons was required to be within 2.5 a of the tt° mass, 
where a is the rms mass resolution, about 5 MeV/c 2 . The 7r° candidates were kinematically 
fit to the 7T° mass to improve momentum resolution; they were required to have a minimum 
momentum of 350 MeV/c. 

To detect the decay Df — > wk + , we reconstructed the to in its dominant decay mode: 
7T + 7r~7r° [0. We normalized to — > r]7i + , i] — > 7r + 7r~7r°, because it has the same final state, 
so the relative reconstruction efficiencies should be near unity and many systematic errors 
cancel in the ratio. We used the CLEO Monte Carlo simulation [IT] to determinine the ratio 
of efficiencies: e(um + )/e(7]7r + ) = 0.91 ± 0.03 (statistical error). The difference from 1.00 is 
primarily due to two kinematic cuts applied to the utt + sample that were not applied to the 
r]-K + sample (described below). 

All requirements were chosen to maximize e/\^N, where the detection efficiency, e, was 
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determined from Monte Carlo, and the background level, N, from the data. The latter was 
done using um + combinations near the mass, but excluding a window around the Df 
mass. 

We began the uj and r\ reconstruction by taking pairs of oppositely charged pions, together 
with a 7T°, and calculating the invariant mass. Three-pion combinations whose invariant 
mass was between 538 and 558 MeV/c 2 (±2cr around the rj mass) were used as rj candidates. 
Combinations with invariant mass between 762 and 802 MeV/c 2 were used as uj candidates; 
this is about a ±0.9 FWHM cut around the uj mass The uj line shape is the convolution of 
its natural width (r = 8.4 MeV/c 2 0) and the detector resolution (<r ^ 8 MeV/c 2 ). 

The rj and uj candidates were combined with a charged pion to form candidates. The 
three charged tracks, two from the rj or uj, along with this "bachelor" pion, were required 
to be consistent with coming from a common vertex. The tracks were refit to pass through 
this vertex, which improves the Df mass resolution by about 4%. 

To take advantage of the hard fragmentation of continuum charm, we required all Df 
candidates to have x > 0.6, where x is the scaled momentum: x = p/p m ax and p ma x = 
i^beam ~ M D + 2 ) 1 ' 2 . This suppresses combinatoric background. A cut on the decay angle of 
the .D+ was also applied. The decay angle (6^) is defined as the angle between the bachelor 
pion in the rest frame, and the momentum in the lab frame. Since the has 
J = 0, the decay angle must have a flat distribution for the signal, while the background 
peaks toward cos 6^ = —1. A cut of cos 6^ > —0.85 was used; this retains 92% of the signal 
and 60% of the background. 

Two kinematic cuts were applied to the ujtt + combinations. First, because the uj is a 
vector particle, it must be produced in the heli city-zero state in the decay — > ujti + . We 
define the helicity angle, a, to be the angle between the normal to the uj decay plane and 
the .D+ direction, both measured in the uj rest frame. This angle must have a distribution 
proportional to cos 2 a. We required | cosa| > 0.45. This cut keeps more than 90% of the 
signal and about 55% of the background. 

Second, the amplitude for the uj decay is maximal at the center of the Dalitz plot. We 
calculated a parameter which is proportional to this decay amplitude; it is simply the cross- 
product of two of the pions' momenta, measured in the uj rest frame. The parameter (R) was 
normalized so that it equals one at the center of the Dalitz plot, and goes to zero at the edge. 
We required R 2 > 0.2; this retains 97% of the signal and about 80% of the background. 

Finally, we sorted the candidates into two categories: "tagged" and "untagged". 
The tagged events are those that are consistent with coming from the decay D* + — > Dfj. 
To tag events, we combined the candidates with photons and calculated the invariant 
mass of each Df'-f combination. To suppress mistags from energy clusters produced by 
hadronic interactions, we required the tagging photon's energy be at least 250 MeV and 
its lateral shape to be consistent with an electromagnetic shower. We calculated the mass 
difference AM 7 = M(D+ 7 ) - M(D+). The D+ is "tagged" if 134 MeV/c 2 < AM 7 < 
154 MeV/c 2 . Events in which no photon meets this criterion are "untagged". The invariant 
mass distribution of the tagged r]7c + combinations is shown in Fig. la. The histogram has 
been fit with a Gaussian for the — ► events and a second-order polynomial for the 
combinatoric background. The mean and width of the Gaussian were fixed to the values 
predicted by Monte Carlo. The fit finds 48. 4^f^ signal events (statistical error only). The 
overlayed functions shown in the figure are the result of a more constrained fit described 
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FIG. 1. Histogram of (a) r]7r + , and (b) um + invariant mass for tagged events. The points with 
error bars are the data; the solid lines are the result of the constrained fit to the data, as described 
in the text. 



below. About 3% of the events contained more than one i]7r + combination which satisfied 
our criteria. The same is true in the ujtt + mode. Since this occurred at the same rate in the 
data and Monte Carlo, and in both the signal and normalizing modes, we accepted these 
double-counts; they have negligible effect on our results. 

A histogram of the invariant mass of the tagged ujtt + combinations is shown in Fig. lb. 
It was fit with the same functions as the i]tt + data, using the same Gaussian parameters, as 
predicted by the Monte Carlo. This fit finds 35.7t];Q;2 — > um + events (statistical error 
only). We consider this to be a significant signal and describe further tests of the data below. 

A number of checks have been performed to help validate this signal. Three-pion com- 
binations were selected in sidebands to the u signal region: 670 MeV/c 2 < M(tt + 7t~7t°) < 
710 MeV/c 2 and 855 MeV/c 2 < M(n+n-n°) < 895 MeV/c 2 . When these are combined with 
a fourth pion, and the ujit + selection criteria applied, no Df signal is seen in either sideband. 
To reproduce the observed um + signal would require a 6 standard deviation fluctuation. 

One can also fit the AM 7 distributions for a signal. Requiring that the four-pion (r)7r + 
or cg>7t + ) mass be between 1943 and 1991 MeV/c 2 and removing the cut on AM 7 , we found 
501 9° r]iT + events and 421^3 um + events, in good agreement with the yields found in the 
previous fits to the r/7r + and um + mass distributions (Fig. 2). In these AM 7 histograms, 
double-counting occurred at a rate of about 10%; this is neglible compared to the statistical 
errors. 

To confirm that these events are in fact — > utt + , rather than some other four-pion 
decay of the D+, we loosened the u mass cut and took all tt + 7t~7t° combinations with masses 
between 650 and 900 MeV/c 2 . These were then combined with a fourth pion; the four-pion 
combinations that passed the tagging criteria (and all other cuts) were kept. Again requiring 
that the four-pion mass be between 1943 and 1991 MeV/c 2 , we made a histogram of the 
three-pion invariant mass. A fit to this histogram yields 44 ± 12 events. However, there are 
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FIG. 2. Histogram of AM 7 for (a) r/TT + events, and (b) lott + events. The points with error 
bars are the data. The solid lines are fits, using a modified Gaussian for the signal, whose shape 
was fixed using Monte Carlo events, and a third-order polynomial for the background. 

also real u's in the ujtt + random combinations under the Df peak. To account for this, we 
performed a sideband subtraction, using upper and lower sidebands in four-pion mass. After 
the subtraction, a fit to the three-pion invariant mass found 32 ± 12 cj's, consistent with our 
previous results. 

We have calculated the invariant mass of the "other" three pion combination in each 
l)7t + candidate event. We define M 3 ' to be the invariant mass of the bachelor tt + with the 
7r~ and 7T° from the u. For the un + events, all of the events in the signal region have 
M 3 ' > 1100 MeV/c 2 . Thus these events are not simply — > i]tt + or — > 07r + (with 
cf) — > 7T + 7r~7r°) events feeding into wk + by combining the pions in the "wrong" order. The 
M 3 ' distribution agrees with the signal Monte Carlo prediction. 

Similarly, we reconstructed events in the signal region as K~tt + 7t + 7t° , as might come from 
D* + decay, by assigning the kaon mass to the negatively-charged track. We found that the 
invariant mass for this alternate particle assignment in every case is more than 2040 MeV/c 2 , 
so these cannot be misreconstructed D* + events. Again, the measured distribution agrees 
with the Monte Carlo prediction. 

The untagged sample of uir + events contains a small excess at the mass (Fig. 3). A 
fit yields 133 ± 57 signal events. Fitting the untagged t]tt + distribution finds 312 ± 31 signal 
events. We included these untagged events in the branching ratio measurement. 

The ratio of reconstruction efficiencies, e(um + )/e(riiT + ), is the same for tagged and un- 
tagged events, so the raw ratio of signal events should also be the same in both samples. For 
the tagged events, we find a ratio of 0.74 ±0.25 ujit + event per r]n + event. For the untagged 
events, the ratio is 0.43 ± 0.19. The two ratios are statistically consistent. 

We also performed a simultaneous fit to the four distributions (r]7i + and um + , tagged 
and untagged), and constrained the ratio of wk + to t]tt + events to be the same for both 
samples. This yielded a ratio of O.SGlg;^; the \ 2 °f the fit to the four histograms was 146.8 
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FIG. 3. Histogram of (a) 7yvr + , and (b) um + invariant mass for untagged events. The points 
with error bars are the data; the solid lines are the result of the constrained fit to the data, as 
described in the text. The fits include a Gaussian of fixed mean and width for D + events near 
1869 MeV/c 2 . In the lower plot, the dashed line shows the background function underneath the 
Df signal; the y-axis scale has been zero-suppressed. 



with 161 degrees of freedom. We used the result of this constrained fit in the branching ratio 
calculation; the fit functions shown in figures 1 and 3 are also the result of this fit. Refitting 
the histograms with the number of — > uj-k + events fixed to be zero yielded a x 2 of 166.9, 
an increase of 20.1. 

Using the ratio of efficiencies determined from Monte Carlo and the r] and uj branching 
fractions to 7r + 7r~7r° H, we determined the branching ratio: 



r(D+ -> ^ 



0.16 ±0.04 ±0.03. 1 



The first error is statistical; the systematic error is dominated by variations in the branching 
ratio caused by varying the cuts used in the analysis. These variations help gauge the 
accuracy of our event simulation. The systematic error also includes contributions from the 
uncertainty in the efficiencies, the branching fractions of the r\ and u, and from variations 
in the result using different fitting functions. 

In order to calculate an absolute branching fraction for — > ujti + , we used the new 
CLEO measurement T(D+ -> r)7r + )/T(D+ -> 0tt+) = 0.47 ± 0.07 §, and the PDG value of 
B(Df -> 07T+) = 0.036 ± 0.009 0. This yields a branching fraction: 

B(D+ -> = (2.7 ± 1.2) x 10~ 3 , (2) 

where all the errors have been added in quadrature. Thus we have observed the decay 
— > co>7r + , which may be the result of annihilation decay, final state interactions, or both. 
We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
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